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ABSTRACT 

An extreme case of electron shock drift acceleration in low Mach number 
collisionless shocks is investigated as a plausible mechanism of initial acceleration 
of relativistic electrons in large-scale shocks in galaxy clusters where upstream 
plasma temperature is of the order of 10 keV and a degree of magnetization is not 
too small. One-dimensional electromagnetic full particle simulations reveal that, 
even though a shock is rather moderate, a part of thermal incoming electrons 
are accelerated and reflected through relativistic shock drift acceleration and 
form a local nonthermal population just upstream of the shock. The accelerated 
electrons can self-generate local coherent waves and further be back-scattered 
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toward the shock by those waves. This may be a scenario for the first stage 
of the electron shock acceleration occurring at the large-scale shocks in galaxy 
clusters such as CIZA J2242. 8+5301 which has well defined radio relics. 

Subject headings: Acceleration of particles — Galaxies: clusters: general — 
Plasmas — Shock waves — Radio continuum: galaxies — Relativistic processes 



Introduction 



In galaxy clusters, the presence of relativistic el ectrons with energies of around GeV 
has been implied by ra. dio synchrotron emissions (e.g.. IWillsonlll970l : iGovoni fc Ferettill2004j : 
Fujita fc Saraziru 120011 ) . Although some possible acceleration mechanisms for those rela- 
tivi stic electrons have been proposed, they are still controversial topics (for a recent review. 



see 



Ferrari et al .1120081 ). The diffuse radio emissions in galaxy clusters are categorized into 



several types such as jet robes from active galactic nuclei, radio halos located at the cen- 
ter of galaxy clusters, and radio relics located at the cluster periphery. In particular, the 
origin of radio relics has been a major mystery. Recent observations suggest that they 
are associated with large-scale shocks which are also thought to be possible sources of the 
ultra- high-energy cosmic rays above ~ 10^*'^ eV ( iKang et al.l 119971 : llnoue et al.l 120071 ). The 
existence of shocks in the radio relic implies that the high -energy electrons are most likely 



prod u ced through the diffu si ve shock acceleratio n process (lEnfilin et al.lll998l : iMiniati et al. 



200ll : iGabici &: Blasil 120031 : iKeshet et al.l |2004| ). In fact, radio shells showing the spec- 
tral softening due to the synchrotron and t he inverse Compton coolings are observed in 
CIZA J2242. 8-1-5301 (|van Weeren et al.l 120101 ). while a temperature jump is found in A3667 
( jFinoguenov et al.ll2010l ). Observations suggest that the temperature is T ^ 9 keV, the 
magnetic-field strength 5^5 /xG, and the Mach number M ^ 4.5 for CIZA J2242. 8+5301, 
and T ^ 2 — 5 keV, -B > 3 yuG, M ^ 2, and the number density is n k, 10~^ — 10"^ for 
A3667. 

According to the standard theory of structure formation, galaxy clusters result from 
mergers of smaller structures. High Mach number shocks (M > 10) are formed in accret- 
ing small halos and around the virialized core of the galaxy cluster. They can accelerate 
cosmic rays to ultra-high energies with harder energy spectrum. On the other hand, low 
Mach number shocks (M < 10) are formed in the virialized core. During the mergers and 
accretions, the most kinetic energy of the accreting rnatter is dissipated by relatively low 
Mach number shocks with M < 4 jMiniati et al.lboool : IRvu et all [20031 iKang et al.ll2007h . 
These shocks with M ^3 — 4 are needed in order to explain spectral indices of the radio 



relic (IGabici fc Blasill2003[ ) . In fact, the shock velocity in the radio relic of the galaxy cluster 
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i s of the order of 10^ km and the correspond ing Mach numbers are about M ^ 2 — 4.5 
( Ivan Weeren et al.ll2010l : iFinoguenov et al.ll2010l ). The shock velocity is similar to what is 
observed for young supernova remnant (SNR) sho cks which generat e galactic cosmic-ray 



electrons (iKoyama et al.lll995l ) as well as nuclei (e.g. lOhira et al.ll201ll ). However, the Mach 
number of a galaxy cluster shock is much smaller than that of a SNR shock because of high 
upstream temperature as well as the so-called magnetization parameter, a = Bf/ATmimeC^ . 
Here, Bi,ni,me and c denote upstream magnetic field, plasma density, electron rest mass, 
and the speed of light, respectively. Then, for the galaxy clusters, T ^ 1 — 10 keV and 
a ~ 0.1 with Bi = 3 fiG and rii = 10^^ cm^^. Such small values of Mach numbers are of the 
same order as those of the shocks observed in the heliosphere, i.e. the terres trial bow shock 
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The diffusive shock acceleration is one of the most plausible a cceleration mechanisms of 



charged particles at a shock front (e.g.. iBlandford fc Eichlerlll978l ). where a resultant energy 



spectrum of the high-energy particles obeys a power law. For SNR shocks, galactic cosmic 
rays are thought to be accelerated through this process. It is not resolved, however, how 
a part of particles originally composing a thermal background plasma is embedded in the 
diffusive shock acceleration which requires pre-existence of nonthermal particles possible to 
flow freely over the shock. This is the so-called injection problem on the diffusive shock 
acceleration. In general, electron injection is thought to be more difficult than ion injec- 
tion because of cross-shock potential. Neverthele s s, relativistic electrons are believed to b e 



accelerated in SNR shocks (IKoyama et al.l Il995l : iBamba et al.l l2003l : iPannuti et al.ll2010l ). 



Although a variety of theoretical as well as simulation studies have revealed that microinsta- 
bilities in the shock transition region play some crucial roles, the electron injection process 
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In the above context, electron injection in low Mach number shocks is inferred to be even 
harder, since the free energy for microinstabilities in the shock transition region should be 
much smaller than that in high Mach number shocks. Althou gh a variety of microinstabilities 
get excited in some heliospheric shocks (e.g., IWu et al.lll984j ). it is known, for instance, that 
the saturation level of the modified two-stream instability which is one of the most plausible 
instabilities in s hocks observed n ear the earth is not so high as that of the Buneman instability 
in SNR shocks (IMatsukiyd |20 lOl ) . Given this perspective, an injection process of electrons in 
the galaxy cluster shocks is veiled in mystery. Indeed, this problem has never been addressed 
so far. Here, we take notice that upstream temperature of the cluster shocks is much higher 
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than those of SNR and heliospheric shocks. As can be seen in the next section, the high 
upstream temperature enables some electrons to be accelerated to relativistic energies via 
the so-called shock drift acceleration (SDA). 

In this paper we investigate a limiting case of the electron SDA and associated micro- 
processes which may be responsible for the injection of electrons into the diffusive shock 
acceleration in the galaxy cluster shocks. In section 2, the SDA process is briefly reviewed 
and extended to an extreme case. Simulation settings and results are shown in section 3. 
Then, discussions and summary are given in section 4. 



2. Shock Drift Acceleration 



2.1. Some Basic Features 



Shock drift acceleration (SDA) is one of the efficient acceleration mechanisms of charged 
particles in the transition region of an oblique coUisionless shock. The motion of non 



relativ i stic electrons in this process is analyzed in detail previously (eg. iLeroy fc Mangenev 



f l984[ ): IWulfll984l ): lKrauss-Varban et al.l (119891 ) : lKrauss-Varban k Wul (119891 ) : lKrauss-Varban fc Burgess 
(jl99ll )). In the so-called normal incidence frame (NIF), where upstream flow direction is par- 
allel to a shock normal, the acceleration occurs while an electron stays in the shock transition 
region and drifts along the shock surface due to finite gradient of magnetic field strength. The 
electron gains energy since the direction of the drift motion is anti-parallel to the motional 
electric field. 

Assuming that the coplanar magnetic field for a one- dimensional shock is in the x — z 
plane and the shock normal is along the leads to 
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-V ■ E 



-X— + VyEo 



(1) 



Here, 7, p, v and e are the Lorentz factor, momentum, velocity and elementary charge, 
E electric field, and the dots denote time derivative, respectively. Furthermore, E = 
{—d(p/dx,Eo,0) has been used, where (p aiid Eq denote shock potential and the motional 
electric field. This results in 



W = 



d 
It 
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Some incoming electrons gain significant energy during the drift motion and are refiected 
backward from the shock. These refiected electrons become a non-thermal component in the 
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upstream plasma frame. The reflection occurs due to the magnetic mirror effect. The de 
Hoffmann- Teller frame (HTF) where upstream flow is along the magnetic field is convenient 
to consider motion of a particle, because the motional electric field disappears in this frame. 
The particle velocity in the HTF (v') and the NIF (v) are related as v' = v + Vht? where 
Vht = (Vi X Bi) X n/(Bi ■ n) is the de Hoffmann- Teller velocity, fi is the unit vector 
normal to the shock, Vi and Bi are upstream flow velocity and the magnetic field in the 
NIF, respectively. In the HTF the right hand side of equation ([2]) becomes zero so that W 
is conserved. Moreover, if the motion of the electron is assumed to be adiabatic, magnetic 
moment, fi' = p']_/2meB\ is also conserved. Here, the prime denotes a quantity measured 
in the HTF. This assumption is valid while the spatial scale of a shock transition region is 
sufficiently larger than the Larmor radius of the electron. From these two restrictions, in 
order for the electron to be mirror reflec ted, the initial pitch angle, a'l = tajLi~^{p[j_/p^^), 
should satisfy the following condition (cf. iFeldman et al.l ( 119831 ) for nonrelativistic version). 



a[ > sin 



B' 



u 



12 



(3) 



Here, B' denotes magnetic field strength, u' = p'/rUeC, A<l>' = e(02 — (f)[)/mec'^, and the 
subscripts 1 and 2 indicate upstream (initial) and downstream quantities, respectively. If 



there exists an overshoot, however. 



i?2 and 



should be defined there. 



The energy gain in the SDA process is obtained by integrating equation in the 
NIF, nieC^Wit) = —e f Eody, as the energy gained by an electron during its drift in the 
transition region along the d irection anti-parallel to the motional electric field (for details to 
Krauss-Varban fc Wul (119891 )). If one move s in the H TF, the energy gain is explained by the 
well known Fermi acceleration mechanism (jWulll984j ). If a particle velocity before reflected is 



written as Vpre, Vpre 



Vpre + Vht- After being reflected, the parallel velocity is reversed so 
that v'post = (~^||,pre' ^^J-.pre)- Then, moving back in the NIF results in Vpost = Vp^st - Vht- 
Now, Vht = (0, 0, f/i tan 6b„) when Vi = Uin and fi is along the x-axis. Hence, the 
acceleration occurs mainly in 2:-directio n. The more Vht inc r eases , the larger energy gain 
of the particle becomes. According to iKrauss-Varban et al.l (119891 ) . the acceleration time 
through this process is typically of the order of an ion gyro period, ~ ^r^' which will also 
be confirmed in our simulation later. For typical intra-cluster magnetic field of fiG, this 
corresponds to about a few 100 seconds, which is of course much shorter than any time 
scales of incoherent loss processes. 

If A$' = is assumed, a[ > a'^ = sin^^ {B[/ B'2y/'^ is the necessary condition for an 
electron to be reflected. When a magnetic overshoot is neglected for simplicity, B2 or the 
critical pitch angle called a loss-cone angle, a'^, can be determined from the Rankine-Hugoniot 
relation by giving a Mach number (Ma), an upstream plasma beta (/3 = 87mi(Te + Ti)/Bf), 
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and a shock angle (8 Bp) wh ich is an angle between the shock normal and upstream magnetic 
field ( Tidman fc Kralllll97ll ). Top two panels in Figure[T]show regions in Ma—Qbu parameter 
space where some electrons on a velocity shell with a radius of 3vte in an upstream plasma 
frame of the NIF can be refiected, where Vte is upstream electron thermal velocity. In the two 
panels, for instance, in the black area some electrons on such a velocity shell corresponding 
to upstream electron beta (/3e = STiriiTe/ Bf) being equal to 0.1 satisfy the condition a'l > a'^, 
where Tj = Tg has been assumed to calculate a'^. Similarly, in the dark (light) gray area some 
electrons can be additionally refiected by assuming f3e = 0.5(1.5). The top (middle) panel 
corresponds to the case with a = 10^^(3"^). The figure confirms that electrons are hardly 
refiected when and/or Qbu become too large. This may be graphically understood as 
follows. Suppose a velocity shell in Vz — fx space indicated as the broken circle in the bottom 
panel. Its center corresponds to the bulk velocity of the upstream plasma in the NIF, [Ui, 0), 
and the radius of the shell is Vr- Let us consider whether some electrons on this velocity shell 
is mirror refiected or not (In this schematic picture relativistic distortion of the shell has been 
neglected.). If Vht is added in Vz, the center of the shell shifts along the thick gray vertical 
arrow. Now a broken line connecting the origin with the new center of the shell is parallel to 
an upstream bulk velocity in the HTF which is along the upstream magnetic field, Bi. Two 
solid lines symmetrical with respect to the broken line denote a loss-cone. Therefore, the 
electrons outside this loss-cone indicated by the black solid arcs can be refiected, otherwise 
transmit. When Ui or a Mach number increases, the center of the shell shifts along the 
broken line and finally a whole shell lies inside the loss-cone. In this case no electrons are 
refiected. Similarly, when G^^ becomes large, the broken line gets more vertical and the 
center of the shell walks away from the origin. Then the whole shell again enters a loss-cone. 
The possible refiection area indicated in the top two panels are defined like this by assuming 
the shell radius Vr = 3vte- The area expands with increasing /3e as expected. On the other 
hand, it does not depend much on a. In particular when a < 10^^, distributions of the 
areas are almost exactly the same as the top panel. It is mainly due to the light speed limit 
of a shell radius that the areas look contracted for a = 3^^ in the middle panel, because 
the increasing a with constant /3e results in increase of the thermal velocity. (The actual 
calculation is performed by taking relativistic effects into account so that the shell is given 
in momentum space with Pr = Spte-)- It is speculated that in extremely high Mach number 
shocks like a SNR shock whose Mach number is typically M4 ~ 10^~^ electron refiection 
hardly occu rs unless there are no preheating mechanisms as discussed by lAmano &: Hoshino 
(120071 . l2010l ). On the other hand, electrons can be relatively easily refiected when Ma < 10, 
although acceleration may be weak because of small Vht- In p ractice, one should note that 
finite shock potential may reduce the refiection rate (jWulll984j ). 
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2.2. An Extreme Case 



As mentioned above, in low Mach number shocks {Ma < 10) electrons are relatively 
easily reflected, although the resultant acceleration is not so efficient in general. However, if 
there are some electrons on the velocity shell with sufficiently large Vr, they may be reflected 
even at large G^n where Vht ~ C){c). The reflected electron energy in such a case would 
become relativistic. 

A condition that an upstream electron on a velocity shell with radius Vr can be reflected 
is vl > (^a/V^t + ^1 ) sina^. When Vht > Ui, this leads to v'^ > Vht\/B[/B^, where 
effects of shock potential has been again neglected. If Vr = rjvte, the following condition is 
finally obtained: 

f/itanGs^ fBi/B2 



V> Vc 



1/2 



(4) 



c V /3eCr/2 

The inequality (jl]) implies that reflected electrons can be present if there are some electrons 
on a velocity shell with Vr > f]cVte- The reflected electrons will have relativistic energies when 
Vht = Ui tan0B„ ~ 0{c). In this limit rjc — )■ r/r, where 



Vr 



B1/B2 



1/2 



a/ Bi/ B2 

Vte/C 



(5) 



In the solar wind or interstellar plasmas rjr is extremely large, since electron temperature is 
~ 100 eV at the highest and a ~ 10~('^~^) (/3e ~ 0(1)). Therefore, satisfaction of the above 
condition for a shock in such an environment as in the earth's bow shock or SNR shocks 
is almost hopeless. However, that may be possible if the electron temperature becomes ~ 
a few— 10 keV or a ~ 10~^^~^^ and /^^ ~ 0.1 — 10 as in some large- scale shocks in galaxy 



clusters flFuiita et all 120071 : iNakar et al.ll2008l : Ivan Weeren et al.ll2010[ ). 



3. ID PIC Simulation 

In this section one-dimensional particle-in-cell (PIC) simulation is performed to repro- 
duce relativistic SDA of electrons discussed above in a self-consistent manner. 

In the simulation a shock is produced by the so-called injection or reflecting wall method. 
An upstream magnetized plasma is continuously injected from the left-hand boundary. The 
plasma is reflected at the right-hand boundary and mixture of the incoming and the reflected 
plasmas results in a downstream medium. The shock is produced at a boundary between 
the upstream and the downstream plasmas. Since the simulation frame is the downstream 
rest frame, the shock propagates in time from right to left. The simulation is done in the 
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NIF so that the injection flow is parallel to the shock normal which is along the x-axis. The 
upstream magnetic field is in the x — z plane. The size of a spatial grid is Ax ~ 0.24ADe 
where Xoe denotes the electron Debye length, and the number of super-particles per cell is 
Np = 200 for both electrons and ions. Time resolution is At = 0.06875ci;~g^ where Upe is the 
electron plasma frequency. Other physical parameters are shown in Table [H For Runs A-C, 
all injection parameters are common except for the shock angle, Qbti- The injection Alfven 
Mach number is M^j„ = 5 which results in ~ 7.8 or Ui ~ 0.061c in the shock frame 
for Run A, for instance, where Ui denotes the upstream bulk velocity. This value roughly 
corresponds to the fast Mach number Mf ^ 2.9 by using a definition of ion acoustic velocity 
as Cg = (3Tj + Tf.)/mi. The magnetization parameter is a = 1/9 and the plasma beta is 
/3 = 3(/3e = (3i) so that the upstream electron temperature is ~ 43 keV. An ion-to-electron 
mass ratio mi/nif. = 1836 is realistic. In the last two columns the maximum value of the 
magnetic field roughly measured in the overshoot relative to the upstream value, Bos/Bi, 
and associated rjc calculated from equation (jl]) are denoted for reference. Note that rjc may 
in fact be larger if non- negligible potential effects are taken into account (jWulll984j ). Most 
efficient electron acceleration is observed in Run A as expected from that [/itanG^n/c is 
the largest, although in all cases reflected electrons are seen. Run D is performed with more 
reahstic parameters where MAin{= 3) and cr(= 1/16) are reduced. The results of Run D are 
essentially the same as those of Run A and discussed in section 4. In the following, results 
of Run A are focused. 



3.1. Run A 



In Figl2] spatio-temporal evolutions of and flelds are shown in the left panels. 
Structure of the shock is more or less time stationary, although it represents weak breathing 
features (jComisel et al.ll201ll ). An identical trajectory of one of the reflected electrons is 
denoted as the black solid lines. Its energy (upper) and momentum (lower) time histories are 
plotted in the right panels. Here, a rate of output data points has been signiflcantly reduced, 
otherwise Px (blue) and Py (green) lines flU in the area because of rapid gyro motions. The 



Table 1: Shock parameters 





Ma 


Mf 


Ui/c 


a 




Te (keV) 


nii 






Bos/Bi 




Run A 


7.8 


2.9 


0.061 


1/9 


1.5 


42.6 


li 


m 


85 


6 


0.99 


Run B 


7.7 


2.9 


0.060 


1/9 


1.5 


42.6 


U 


^36 


80 


6 


0.48 


Run C 


7.4 


2.8 


0.058 


1/9 


1.5 


42.6 


U 


^36 


60 


6 


0.14 


Run D 


5.0 


1.9 


0.029 


1/16 


1.5 


24.0 


U 


^36 


87 


3.5 


1.4 
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electron gains energy mainly through pz which is almost parallel to the magnetic field during 
its stay in the shock transition region. The time the electron stays in the transition region 
during the reflection proces s is ~ VLT^. These are typical features of SDA discussed by 



Krauss-Varban &: Wul (Il989l ) and iKrauss-Varban et al.l (119891 ) , while energy of the reflected 



electron here becomes relativistic. 

The flrst three panels from the top in Fig|3]show ion v^—x, electron p^. — x, and electron 
Pz — X phase spaces at ojpet = 64350 (fijt ~ 11.7) which is indicated as the dashed lines in 
FigJ2J It is clear in the third panel that some electrons are reflected basically along the 
magnetic fleld and have relativistic energies. The fourth panel represents p± — p\\ phase 
space of the electrons surrounded by the black square in the third panel. Their energy 
distribution function is plotted as the black line in the bottom panel. The reflected electrons 
show a ring-beam feature leading to the non-thermal part of the distribution function. A 
fraction of the bulk energy density of the incoming ions carried by the reflected electrons, e, 
is estimated as e ~ 0.014. The gray line in the bottom panel is a downstream distribution 
function corresponding to the region surrounded by the gray square in the third panel. The 
high energy part (7e — 1 > 1) looks nonthermal. Its origin is the electrons also having large 
negative Pz around x/pi ~ 11 and 13 in the third panel. They are produced in the second 
and the third magnetic overshoots at x/pi ~ 12 and 14 seen in Figl2l By comparing the 
upstream and the downstream energy distribution functions, the highest energy electron is 
upstream. This implies that in the simulation the downstream region is still very limited to 
see an equilibrium state and that is clear also from the top three panels. 

It is often thought that reflected electrons produce a loss-cone distribution because 



of their adiabatic behaviors (ILobzin et al.ll2005l ). However, no clear loss-cone is seen in 



the fourth panel. This is probably due to the effects of small-scale waves generated in 
the transition region. Although we avoid getting involved here with detailed analysis of this 
problem, we just show evidence of the presence of such small-scale waves. Fig 0] shows By fleld 
fluctuations in a transition region (upper panel) and its Fourier spectrum at Upet = 44220 
(lower panel) indicated by the solid line in the upper panel. Tiny streaks in the transition 
region are visible in the upper panel. They appear in the lower panel as spectral peaks 
between 0.1 < kX^ < 1 where Ag = c/upe denotes the electron inertial length. Corresponding 
structures are also seen in ion Vy — x phase space indicating that reflected ions destabilize 
some kind of microinst ability (although not shown). 
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3.2. Injection of Reflected Electrons into Further Acceleration Process 

Here, some possible behaviors of the accelerated electrons after being reflected are dis- 
cussed. Since the reflected electrons form the non-equilibrium distribution function, some 
local instabilities may be driven. However, they may not have been properly treated in 
the previous section, because there wave vectors only along the x-axis are allowed. In the 
following possible instabilities and associated wave-particle interactions are discussed by 
performing further one-dimensional simulations for various wave propagation angles with 
periodic boundary conditions. A plasma is assumed to be composed of three components, 
i.e., incoming electrons and ions, and reflected electrons. The reflected electrons are assumed 
to form a ring-beam distribution with a bulk momentum {j)\\o/mf.c,p ±^0/171^0) = (2, 1) and a 
relative density nb/ui = 0.1, where rii denotes a density of the incoming ions which are at 
rest in average. Hence, the simulation is done in the upstream plasma frame of the reference. 
A density and bulk momentum of the incoming electrons are decided to satisfy charge and 
current neutral conditions. For all three components, initial temperatures are equal and 
isotropic in their proper frames, T/mgC^ = Pe<^/2, where Pe and a are chosen to be same as 
the values of the upstream plasma in the previous section so that Pe = 1-5 and cr = 1/9. 
The system size is L/Xe = 819.2, the number of spatial grids = 8192 corresponding to 
Ax = O.lAg ~ 0.35\r)e, the number of super-particles per cell Np = 400, and a time step 
ujpeAt = 0.1, respectively. With these parameters flxed, four different runs are performed 
with various wave propagation angles {Osk = 0°, 30°, 60°, and 80°) with respect to the mag- 
netic fleld which is in the x — z plane. Note that Obr is independent of the shock angle, 

©Bn- 

Before discussing simulation results, we flrst briefly note how the number of super- 
particles per cell, Np, is crucial in reproducing a long time evolution of the system. FigJS] 
represents evolutions of spatially averaged magnetic fleld energies as a function of time in 
0Bk = 60° case for three different values of Np] thin black lines correspond to Np = 4, thin 
gray lines to Np = 40, and thick black lines to Np = 400, respectively. The solid and the 
broken lines indicate By and components. Because of the large noise level, Np = 4 causes 
totally different results from others. For Np = 40, early development of an instability appears 
to be calculated properly. However, nonlinear evolution of the system {upet > 1000) may 
not be well described. In contrast to the case with Np = 400, the Bz component settles in 
a constant value after Up^t = 1600 and this results in the same level of magnitude as the 
By component at last. For such an oblique propagation angle, B^ can easily couple with 
electrostatic fluctuations. The lower limit of B^ fleld energy in this nonlinear stage may be 
an influence of electrostatic noise. We conflrmed that the results are qualitatively common 
for Np > 100 at least up to the time to be discussed here. In the following, results for 
Np = 400, which we believe to be reliable, will be discussed. 
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Figj6] represents time histories of field energies (top) and effective electron temperature 
anisotropy (second), u — k spectrum of field for < Up^t < 204.8 (third), and electron 
distribution in — py phase space at UpJ: = 3200, respectively, for Osk = 0°. Here, the 
effective electron temperature anisotropy is defined as Te\\/Te± = T.jUjp'^^^/T.jnjp'j^, where 
the summation is taken over all the incoming and the reflected electrons. The rapid growth 
of Ex field energy is essentially due to the beam (two-stream) instability in which the electron 
beam destabilizes mainly Langmuir waves and their higher harmonics as shown in the third 
panel. Although the electron distribution function forms a plateau in p\\ in the end of the 
run as seen in the bottom panel, the temperature anisotropy still persists (second panel). 
Because of the small electromagnetic field energies throughout the run {Ey and E^ field 
energies are in the noise level.), the above process is essentially electrostatic. Electrostatic 
wave activities are again dominant for Osk = 30°, although their intensity is much less than 
the case with Osk = 0° (not shown). 

In contrast, electromagnetic fluctuations become predominant for Osk ^ 60°. The simi- 
lar plot as Figj6]for 9Bk = 60° is shown in FigJTJ although the third panel is u — k spectrum 
of By field for < Upet < 3276.8. This electromagnetic instability is essentially nonresonant 
type which is confirmed from that natural modes of the beam electrons in the u — k spectrum 
(along u ~ kvxo {vxo is a projection of a parallel bulk velocity into the x-direction)) have no 
significant intensiti es in the thi rd panel. They are probably related with one of the oblique 



modes discussed by iBretl ( 120091 ). although detailed linear analysis of kinetic relativistic ring- 
beam instabilities in a magnetized plasma should be reported elsewhere. The nonresonant 
instabilities more efficiently relax electron temperature anisotropy (second panel) than the 
electrostatic beam instability. This is also confirmed as the widely scattered ring-beam elec- 
trons in p± — p\\ space in the bottom panel. Especially, some of the ring-beam electrons 
have negative p\\. This means that some of the reflected electrons are back scattered by the 
self-generated waves. (More precisely, electrons having v\\ cos Qbu < ^sh can be regarded as 
being back scattered {vgh denotes the shock velocity), although we have not specified here a 
value of Qbu-) Some typical trajectories of the back scattered electrons are shown as solid 
lines in the top-left panel of FigJH] in which the background gray scale denotes amplitude of 
the magnetic fluctuation |B — Bo|, where Bq is the ambient magnetic field. While all the 
electrons initially propagate in the positive x-direction, they finally have negative velocities 
in X. Time evolution of the pitch angle cosine of the electron corresponding to the black solid 
line is plotted in the top-right panel and its trajectory in p±_ — p\\ space is indicated in the 
bottom panel. Major changes in the pitch angle cosine occur in two bounded time domains 
where the electron encounters large amplitude wave packets {cOpet ~ 500 and ~ 1800). In 
most of the remaining time the electron propagates in one-direction, although rapid changes 
in its pitch angle exist. These features are common for all other trajectories in the top-left 
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panel. 

For Osk = 80°, a nonresonant instability is again dominant (not shown). In such a 
large dsk generated waves are almost nonpropagating pure growing modes which are basic 
features of the Weibel instability. However, growth time is much longer and relaxation of 
the effective electron temperature anisotropy is less efficient than Osk = 60° case. 



4. Summary and Discussions 

It was shown by using a one-dimensional electromagnetic full particle simulation that 
relativistic SDA of electrons is feasible even though a quasi-perpendicular shock is rather 
moderate when the upstream electron temperature becomes of the order of 10 keV or the 
magnetization parameter a is not too small (equation (5)). Such a condition may be realized 
in large-s cale shocks of gal a xy clu sters, cosmic-ray-modified subshocks of SNRs, etc. For 
instance, Ivan Weeren et al.l ( l2010l ) showed an evidence of strong acceleration of relativistic 
electrons forming radio relics in the merging galaxy cluster CIZA J2242. 8-1-5301, where an 
average temperature of the intra-cluster medium is estimated as ~ 9 keV. Their polarization 
analysis indicates that the observed shock in the radio relic is quasi-perpendicular. 

The maximum energy of the reflected electrons is affected by the de Hoffmann- Teller 
velocity, Vht = UitanQBn, which is a strong function of a shock angle. The Lorentz 
transformations derive energy of a reflected electron as jref = l^'jiTi^ + (^ht ~ '^f^HTVz/c), 
where 7 and Vz denote the electron's initial NIF Lorentz factor and velocity component in 
z-direction in the shock frame, and Tht = (1 ~ /^ht)"^^^; (^ht = Vht/c, respectively. When 
the nonrelativistic limit with large 9b„,(~ 90°) is considered, the above expression reduces 



to equation (16) in iKrauss- Varban &: Wd (119891 ). Observed values of the maximum Lorentz 



factors of upstream electrons for Runs A-C are ~10.0, ~4.1, and ~2.7, respectively. They 
roughly coincide with the values derived from the above estimate, ~9.6, ~3.6, and ~2.3, 
where 7 ~ 1.9 (at a maximum) and Vz/c ~ —0.8 are again from the simulation. Possible 
reasons for the underestimate may be neglecting a temporal variation of the shock velocity 
and/or nonadiabatic features based on the small scale waves in the transition region. 

The parameter dependence of the maximum attainable energy through the SDA is 
shown in FigJHl Energies of the reflected particles along the edges of the possible reflection 
areas in the top two panels in FigJT] are plotted. Each line color corresponds to different Pe 
(black: 0.1, dark gray: 0.5, light gray: 1.5). The upper panels show energies of the reflected 
electrons normalized to upstream electron bulk flow energies. The dashed lines denote bulk 
energy of the ion flow. The normalized energies indicate apparent M7^ dependence for wide 
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range of Ma- This is consistent with the resuh that the obtained energies do not basically 
depend on Ma as confirmed in the lower panels. This is due to the fact that Vht is roughly 
constant along the edges of the possible reflection area. However, the maximum energy 
strongly depends on a implying that the magnetic field strength is crucial. Indeed, it should 
be noted that the vertical axes in the lower right panel is three orders of magnitude larger 
than those in the lower left panel, while a ratio of a between these two panels is also ~ 10^. 

The reflected electrons are accelerated mainly in 2;-direction which is almost along the 
magnetic field, forming the nonthermal part of the energy distribution function upstream of 
the shock. The efficiency of the energy transfer to the reflected electrons is highest in Run 
A. If an injection rate, e, is defined as the fraction of energy density carried by the refiected 
electrons with respect to the bulk energy density of the incoming ions, e ~ 1.4% in this 
case. The corresponding relative number density of the refiected electrons, or a refiection 



ratio, is 2.6%. These values are to be compared with observations. iFinoguenov et al.l (120101 ) 
estimated the total energy of radio-emitting electrons in the radio relic in A3667 is 0.15% of 
the shock kinetic energy, where they adopt the field strength B ^ 3 fxG to explain dim inverse 
Compton emission in the X-ray band. Another observational estimation of the injection 
efficiency is obtained assuming the equipartition between the energies of the magnetic field 



and th e radio emitting electrons. Using the observed quantities given in IFinoguenov et al. 



(I2OIOI ). we derive 4% of the shock kinetic energy goes into that of radio emitting electrons. 
The injection efficiency of a few percent obtained by our simulation result is comparable 
to the observational estimation (0.1-4%), which, however, contains various uncertainties. 
It should be noted that the result here is originated from the one-dimensional simulations 
only for the particular parameter set, and that the injection efficiency depends on various 
parameters such as the electron temperature, the shock normal angle, etc. To estimate the 
injection efficiency more accurately, we should perform detailed survey on parameters as well 
as two or three dimensional simulations which will be addressed in future works. 

The accelerated relativistic refiected electrons form a non-equilibrium local ring-beam 
like distribution function upstream of the shock. This is possible to generate large ampli- 
tude waves through a variety of microinstabilities. Local simulations with periodic boundary 
conditions reveal that the electrostatic electron beam instability is rapidly destabilized ba- 
sically along the magnetic field. This instability quickly saturates by forming a plateau of 
the distribution function so that the effective temperature anisotropy {Te\\/Te± > 1) remains. 
Afterwards, the oblique nonresonant instability grows slowly leading to efficient relaxation of 
the temperature anisotropy. In this process some of the refiected electrons are scattered back 
toward the shock by self-generated coherent wave packets. A sequence of the above process, 
i.e., the refiection of some incoming electrons through the SDA process followed by excitation 
of the microinstabilities and the backscattering of a part of the refiected electrons by the 
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self-generated waves, may be the first step of injection into the diffusive shock acceleration 
process. 

In Runs A, B, and C, Ma and Tg may be somewhat higher than those of typical large- 
scale shocks in galaxy clusters. Therefore, an additional Run D is performed to confirm 
that the similar process can work also for more realistic parameters. Here, Ma ~ 5.0(t/i ~ 
0.029c), Tg ~ 24.0 keV, and a = 0.0625 as shown in Table [TJ Results are qualitatively similar 
to Run A. The relativistic SDA works and associated ring-beam electrons are produced 
upstream. However, the maximum Lorentz factor of the reflected electrons, 7^^^ ~ 4, is not 
as large as in Run A, since Vht and the initial 7 are a little less because of the small Ma and 
Te here (Fig JTOj) . In addition, a reflection rate of the incoming electrons becomes smaller, 
because rjc ~ 1.4 which is larger than rjc ~ 0.99 in Run A. Nevertheless, the relativistic 
SDA can work in such a parameter regime too. In practice the process may occur even for 
more moderate and lower temperature sho cks if a halo electron component is pres ent in an 



intra-cluster medium as in the solar wind (jStverak et al.ll2009l : iLouarn et al.ll2009l ). 



In the present paper effects of higher spatial dimensions have been excluded. We expect 
that the processes discussed separately in sections 3.1 and 3.2 are simultaneously observed in 
2D or 3D systems. However, details of the competing processes among them have been un- 
known. For example, the electron beam instability is accompanied by pitch ang le scattering 



through multidimensional wave-particle interactions in its nonlinear stage (e.g. iPavan et al. 



20091 ). If this process proceeds more rapidly than the growth of the oblique nonresonant 
instability, the relaxation of the effective temperature anisotropy and associated backscat- 
tering of a part of the reflected electrons might occur through this. It is also curious to know 
some other effects like a rippling and upstream turbulence. All these are the future issues. 
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Fig. 1. — Feasibility of SDA. Top two panels show Ma — ©Bn parameter spaces where 
adiabatic electron reflection is possible for three different upstream electron beta values with 
a = 10~^ (top panel) and cr = (middle panel). The bottom panel shows a schematic of 
a loss-cone in Vz — Vx space. 
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Fig. 2. — Left: Spatio-temporal evolutions of color coded (upper) and (lower) field 
components and a trajectory of a typical reflected electron. Right: Energy (upper) and 
momentum (lower) histories of the traced electron. 
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Fig. 3. — Phase space densities at Upet = 64350. Top three panels represent ion — x, 
electron — x and Pz — x phase spaces. The fourth and the fifth panels show p± — p\\ 
phase space and a corresponding energy distribution functions of electrons surrounded by 
the squares in the third panel. 
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Fig. 4. — Small scale fluctuations of By in the shock transition region (upper panel) and its 
Fourier spectrum at cupet = 44220 (lower panel). 
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Fig. 5. — Time histories of By (solid lines) and Bz (broken lines) field energies for different 
numbers of super-particles per cell, Np, in case of Osk = 60°. The black thin lines, the gray 
thin lines, and the black thick lines denote Np = 4, 40, and 400 cases, respectively. 
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Fig. 6. — Time histories of field energies (top) and effective electron temperature anisotropy 
(second), u — k Fourier spectrum of E^^ field for < Upet < 204.8 (third), and a final electron 
distribution in p± — py space (bottom) for Osk = 0°. 
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Fig. 7. — Time histories of field energies (top) and effective electron temperature anisotropy 
(second), u — k Fourier spectrum of By field for < UpJ, < 3276.8 (third), and a final electron 
distribution in — P\\ space (bottom) for 6Bk = 60°. 
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Fig. 8. — Trajectories of some electrons back scattered by self-generated waves with the 
background gray scale of amplitudes of magnetic fluctuations (top left). An evolution of a 
pitch angle cosine (top right) and a trajectory in space (bottom) of the black particle 

in the top left panel are also shown. 
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Fig. 9. — Parameter dependence of the maximum energy through the SDA. The upper panels 
show the maximum attainable energy of the reflected electrons along the upper edges of the 
possible reflection areas indicated in FigJT] normalized to bulk flow energy of the upstream 
electrons. The dashed lines denote bulk energy of the upstream ion flow. The lower panels 
represent the corresponding Lorentz factor (left axes) and energy in keV (right axes). The 
line colors indicate different /3e (black: 0.1, dark gray: 0.5, light gray: 1.5). The left and the 
right panels correspond to a = 10"^ and 3~^, respectively. 
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Fig. 10. — Electron phase space density in — x at Upet = 85200 for Run D. 



